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Abstract 

We determine the energy eigenvalues and eigenfunctions of the harmonic osciUator where the 
coordinates and momenta are assumed to obey the modified commutation relations [xi,pj] = ih[{l+ 
j3p^)5ij + (3'piPj]- These commutation relations are motivated by the fact they lead to the minimal 
length uncertainty relations which appear in perturbative string theory. Our solutions illustrate 
how certain features of string theory may manifest themselves in simple quantum mechanical 
systems through the modification of the canonical commutation relations. We discuss whether 
such effects are observable in precision measurements on electrons trapped in strong magnetic 
fields. 
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I. INTRODUCTION 



In this paper, we derive the exact solution to the Schrodinger equation for the harmonic 
oscillator when the commutation relation between the position and momentum is modified 
from the canonical one to 

[x,p\ = in{i + (3p') . (1) 

This commutation relation leads to the uncertainty relation 

which implies the existence of a minimal length 

Aa;mm = Ti^fp , (3) 

below which the uncertainty in position, Aa;, cannot be reduced |T], 0. 

The minimal length uncertainty relation, Eq. @, has appeared in the context of per- 
turbative string theory [Q where it is a consequence of the fact that strings cannot probe 
distances below the string scale h^/^. Though the modified commutation relation, Eq. (P, 
has not been so far derived directly from string theory, the fact that it implies Eq. 
suggests that it is one possible way in which certain features of string theory may manifest 
themselves in low energy quantum mechanical systems. 

It should be noted, however, that Eq. (|^) is not an ubiquitous prediction of string the- 
ory. Indeed, both in the realms of perturbative and non-perturbative string theory (where 
distances shorter than the string scale can be probed by D-branes Q), another type of 
uncertainty relation involving both spatial and time coordinates has been found to hold 
1^. The distinction (and relation) between the minimal length uncertainty relation and the 
space-time uncertainty relation has been elucidated by Yoneya [Q. 

Nontheless, Eq. @) does embody an intriguing UV/IR relation : when Ap is large (UV), 
Ax is proportional to Ap and therefore is also large (IR). This type of UV/IR relation has 
appeared in several other contexts: the AdS/CFT correspondence non- commutative 
field theory and more recently in attempts at understanding quantum gravity in asymp- 
totically de Sitter spaces |P, |T0| . 

It has been argued that the kind of UV/IR "bootstrap" described by Eq. (^ is necessary 
to understand observable implications of short distance physics on infiationary cosmology 
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TT|. Likewise, Banks has argued that some kind of UV/IR relation should be at the core of 



the cosmological constant problem [|T^ as well as its relation to the problem of supersym- 
metry breaking. Therefore, both Eq. (0) and the underlying Eq. (|1|) are well motivated by 
a variety of applications, including the cosmological constant problem which we will discuss 
is a subsequent paper 



Furthermore, the UV/IR relation represented by Eq. (j^) suggests that certain "stringy" 
short distance (UV) effects may manifest themselves at longer distances (IR), lending addi- 
tional justification to our analysis of the non-relativistic harmonic oscillator. 

The problem of solving for the energy eigenvalues and eigenstates of the harmonic oscil- 
lator with the minimal length uncertainty relation has been studied previously by Kempf et 
al. in Refs. [jl|, |^. However, the exact result had been obtained only for the 1-dimensional 
case. We present here the exact solution for the general D-dimensional isotropic harmonic 
oscillator. 



II. THE HARMONIC OSCILLATOR IN 1 DIMENSION 



We represent the position and momentum operators obeying Eq. (|lD in momentum space 



by 



X = ih 
p = p . 



(4) 



The choice of the constant 7 determines the weight function in the definition of the inner 
product: 



dp 



;i + Pp^ 



where 



a 



1 



(5) 



(6) 



This definition ensures the hermiticity of x. In the following, we will keep 7 arbitrary, though 
eventually we will find that the energy levels in fact do not depend on its value. 
The Schrodinger equation for the harmonic oscillator with Hamiltonian 



n = —uuj X H p 



(7) 
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is given by 



2E 



A change of variable from p to 



l=tan'\^p) 



VP 



maps the region —00 < p < 00 to 



(8) 



(9) 



TT TT 

<P< 



2^ 



and casts Eq. (|^) into the form: 



Defining dimensionless parameters by 



P 



1 



2E 



we obtain 



where we use the shorthand notations 



^(0 = 



(10) 



-{^ - (1 + j ) tan^ ^/^P + {^ + nv) = • (11) 



(12) 



(13) 



c = cos , s = sin . 



(14) 



Let ^1/(0 = c^^^ /(■s), where A is a constant to be determined. Then the equation for /(s) is 



[l-s')f"-{2\ + l)sf + 



^-A} + {A(A-l)-lj£, 



/ = o. 



(15) 



The variable is now — 1 < s < 1. Note that 5 = 7//? is eliminated from the equation. We 
fix A by requiring the coefficient of the tangent squared term to vanish: 



A(A-l) 



1 



. 



(16) 
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The wave function should be non-singular at c = 0, which implies 



^ 1 1 1 
'=2 + V4 + ;^ 



(17) 



This simplifies Eq. ([151 ) to 

(l-.2)r-(2A + l)./'+(-^-A)/ = 0. (18) 

Similarly, f{s) should be non-singular at s = ±1. Thus we require a polynomial solution to 
Eq. (|TBp. This requirement imposes the following condition on the coefficient of /: 

e 



- X = n(n + 2X' 



where n is a non-negative integer |T5[. Eq. ([18|) becomes 

(1 - s^)f" - (2A + l)s/' + n(n + 2A)/ = 
the solution of which is given by the Gegenbauer polynomial: 

fis) = c::{s) . 

The energy eigenvalues follow from the condition Eq. (|19]): 

'n^ + (2n + l)A 



K 



[n^ + n + -] + {2n + 1) \ 1 + 



4 ' 



or more explicitly, 

E„ = hui 



1 



This result agrees with Kempf pi]. The normalized energy eigenf unctions are: 



n\ (n + A) J/3 , , . , 
27rr(n + 2A) ^ ' 



where 



c = cos 



l+/3p2 

l3p 



(19) 



(20) 



(21) 



(22) 



(23) 



(24) 



III. THE HARMONIC OSCILLATOR IN L> DIMENSIONS 



In more than 1 dimension, the modified commutation relation can be generahzed to the 
tensorial form: 

[xi,pj] = ih{Sij + Pp'^Sij + P'piPj) . (26) 
If the components of the momentum pi are assumed to commute with each other, 

[P^,P,]=0, (27) 

then the commutation relations among the coordinates Xj are almost uniquely determined 
by the Jacobi Identity (up to possible extensions) as [||, 0] 



These operators are realized in momentum space as 



d d 

dpi dpj 



Xi = ih 

Pi = Pi. (29) 



The arbitrary constant 7 in the representation of Xi does not affect the commutation relations 
among the Xj's. Again, its choice determines the weight function in the definition of the 
inner product: 

where 

D - r 



Note that when P' = 0, this expression reduces to that of the ID case, Eq. (|^). Ref. uses 

7 = /3 + /3'(^), (32) 

in which case the weight function is constant. We will keep 7 arbitrary in our calculations. 
As in the ID case, we will find that the energy eigenvalues do not depend on 7. 
Since the D-dimensional harmonic oscillator 
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is rotationally symmetric, we can assume that the momentum space energy eigenfunctions 
expressed in terms of the radial momenta can be written as a product of spherical harmonics 
and a radial wave function: 



(34) 



In 2D and 3i5, we have 



*2(p) 

^3(P) 



1 

— 



-im(f> 



R{p) 



R{p). 



(35) 



Eq. (^4D allows us to make the replacements 



^d^__^ D-i d 

~[ dpi dp^ p dp p2 
^ d d 
~{ opi dp 



(36) 



where |]T4| 

L2 + -2) , i = 0,1,2,- ■■ . (37) 

{i = \m\ for D = 2.) 

We therefore find, not unexpectedly, that the Schrodinger equation for the D-dimensional 
harmonic oscillator can be reduced to the 1-dimensional problem for the radial wave function 
R{p). The equation for R{p) is: 



[I + {(3 + /?')/] 1 1 +l^ + [iD-l)P + 2^]p\l[l + i(3 + (3')p'] I 



L2 

p2 



+ ( - 2pL^ ) + { 7 ( /5D + /3' + 7 ) - P^L^ } p^ 



1 9 E 

R{p) + --p'R{p) = —R{p). 

(38) 



Introducing the change of variable 



maps the region < p < oo to 



< p < 



TT 



2^pTW 

and renders the equation into the following form: 



(39) 



(40) 



^ d^ iD-l)J(3 + (3' [(D-l)/3 + 27] f- — -\d 



dp"^ 



tanJp + p'p 



dp 



tan^ ^fi + P'p 



ip + P') 



tan^ JP + P'p 



+ 



1 



tan' 



y/P + P'p 



nhuj{p + p') 

Defining the dimensionless parameters 



p 



we obtain 

1 d^R 



, K = J phuiP + P') , rj 



P 



+ 



+ 



{D-l)- + {{D-l)ri + 2S}- 
s c 



P + P 



IdR 



7' ^ 



7 



P + P' 



1 1 s' 



1^ - {2r) - 1)L' + 6D^ - - + [6{D - l)7j + 5(1 + 5) - tj'L' - -j} 



(41) 



= , 
(43) 



where we use the shorthand notation 



c = cos , s — sin , 



(44) 



as before. Let R{^) — c^'^^f{s). The equation for f(s) is 



{2X + l + {D-l){l-r])}s + 



D - 1 



/ = o. 



(45) 



Note that as in the ID case, S is ehminated. We choose A to cancel the tangent squared 
term: 



X'-X[l + iD-l)ri]-ri'L'-- = 

AC 



Taking the positive root, we obtain 



2 V 4 ^4 



The equation for f{s) then simphfies to 



(46) 



(47) 



(1 - s')f 



+ 



{2A+l + (L'-l)(l-7^)}s + 



D-1 



f 
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+ 



|^_Ai^-(2,-l)L^}- 



L2" 
T2" 



/ = o. 



(48) 



Next, let /(s) = s^g{s). This substitution eliminates the centrifugal barrier term and gives 
the equation for g{s): 



+ 
+ 



-{2X + 2i+l + {D-l){l-r])}s 



21 + D -I 



- 2t]L'^ - (2£ + D)X + i{{D - l)r] - 1} 



9 = 



Another change of variable 

z = 2s^ -1 

maps the interval 0<s<lto— l<z<l and leads to 
+ 



(49) 



(50) 



(b-a) - {a + b + 2)z 



dz 



+ - 



27]L^ - {21 + D)\ + i{{D - 1)7] - 1} 



9 = 



(51) 



where 



a = A - 
2 



1 + (D - 1)7] 

2 

-e-1 . 



^ {i + {D-i)vr , 1 

4 + 



(52) 



Once again, we need a polynomial solution to 9{z) to keep the wave function regular at 
z = ±1. The condition one must impose is 



— - 27]L^ - {2i + D)\ + i{{D - 1)7] - 1} 



7i'{7i' + a + 6 + 1) , 



where 71' is a non-negative integer [jTSl. This casts Eq. (^) into the form 



2,d 9 



{h- a) - {a + h + 2)z ^ + n'(n' + a + 6+l)5( = 0. 



the solution of which is given by the Jacobi polynomial: 

9{z) = Pi^''\z) . 

The energy eigenvalues are given by 



(53) 



(54) 



(55) 



D 



{2n +i) + ^ A + 4n''' + 2n{l - 7]){D - 1) + 47i'i + 27]L^ + {7]{D - 1) - l}^ 
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2 ( n + ^ ) A + + n{l - r]){D - 1) + {2r] - l)U 



2 ( n + - \ ^ ^ ^ ^ + 772^2 + — 



+ 



+ + — r/ + —(1 - i]) + {2'q - 1)L'^ 



(56) 



where in the second hne we have made the replacement 2n' + i = n. The final expression is: 



E„ 



huj 



2)\1 



n 



1 + < /32L2 + ^ , ) jl^h^OJ^ 



D 



D2 



(57) 



This equation represents the main result of our paper. The ID result can be reproduced from 
this expression by setting = 1, L2 = 0, and /5' = 0. The normalized energy eigenf unctions 



are 



Rni{p) 



\ 



2{2n' + a + b+l)n'\ V{n' + a + b + l) 



T{n' + a + l)V{n +h + l) 
where n' = {n — i)/2, and 



{P + PY/'c^^'/Pi?''\z) , (58) 



c = cos J P + P' p 



I \/P + (3'p 

v/i + (/? + W 
{(3 + (3')p' - 1 



z = 2s^-l 



{(3 + (3')p' + 1 



(59) 



IV. COMPARISON WITH PREVIOUS RESULTS 

Kempf [0 has calculated the energy eigenvalues of the 2D and 3D harmonic oscillators 
to linear order in f] and f3'. From the exact expression Eq. (0), we can easily identify the 
terms to that order to be: 



, (60) 



For ease of comparison, we have introduced the notation 



(61) 



10 



For the D = 2 and D = 3 cases, the exphcit expressions are 



E2D ~ Tioo 



[n 



+ 1) + ^ {(A;^ + k'^n + 1)2 + - k''){f + 1) + k''} 



2 

3\ 1 



n + ^)+^Uk' + k")[n+^) +{e~k 



/2x 



We define the parameter s by 



s = n' + 1 = — h 1 , 



which takes values from 1 to [{n + 2)/2] for fixed n. Then 



E. 



2D 



E- 



3D 



huj 



(62) 



(63) 



{n + l) + k^{n^ + 3n + 3) - k' i^n + -j - 2 (P - k' ) s{n + 2- s) 
n + ^^+e (n^ + 4n + - k'^ + ^) - i^'' - k'^) s(2n + h- 2s) 



(64) 



which agree with Ref. [Q]. See also Ref. ||T6| 



V. ELECTRONS IN A PENNING TRAP 



The TT? dependence of the harmonic oscillator energy levels, Eq. (pT]), would be an un- 
mistakeable signature of the modified commutation relations Eq. (|lD or (pB]). Even if (3 and 
j3' were small, the deviation from the usual n dependence should be manifest for sufficiently 
large n. 

The cyclotron motion of an electron in a Penning trap [|r^ is effectively a one-dimensional 
harmonic oscillator. By looking for shifts in its energy levels, we may be able to place a 
constraint on /3. To leading order in j3 and n, the deviation of the harmonic oscillator energy 
levels from the canonical hu{n + |) is given by 



AEn ( [3mhuj\ 2 

hu ~ \ 2 ; 



(65) 



which grows quickly with n. Note that the combination of parameters that is constained by 
a measurement of AEn/huj is jSmhuj. The cyclotron frequency of an electron trapped in a 
magnetic field of strength B is (in SI units) 

eB 



OJr. 



(66) 
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1-53 i„2 _2 /„2 



Therefore, 

ruehuc = {eK)B = (1.7 x 10"^^ kg^ ■ mVsVT) B , (67) 

which is independent of the electron mass mg. For a trapping field strength of 5 = 6 T, we 
obtain 

rriehuj^ = ehB = (1.0 x 10"^^ kg^ ■ mVs^) • (68) 

Even though we anticipate that measuring the energy levels accurately for very large n would 
be difficult, let us assume for the sake of argument that deviations as large as would be 
detectable. Then, the absence of such a deviation for the n-th energy level would imply 



(^) 



< 1 , (69) 



or 



KenB / 



This translates into 



hJf3 < 



/2h\ 1 _ (1.5 X 10"^ m 
ClS db limit for the minimal length, Eq. (0), and 



- = — , (71) 

n'' n 



1 

> 



'ehB" 



= (7.1 X 10"^^ kg ■ m/s) n = (13 eV/c) n , (72) 



for the string momentum scale. 

Aside from the technical question of whether one can measure the energy levels accurately 
for large n, we must require that n satisfy 

§ « 1 , 73 

for our electron to stay non-relativistic. This leads to the constraint 

This condition also maintains the radius of the electron's cyclotron motion to be well within 
the geometry of the Penning trap. Therefore, the maximum value of n that can be used to 
constrain j3 would be n ~ 10^ if we allow for a 10% relativistic correction. So the best limit 
on P that can be achieved will be 

hJp < 10-^^ m , 4^ > 1 GeV/c . (75) 

VP 

Obtaining a better limit would be difficult since improving the limit on y/]3 by an order of 
magnitude requires the improvement of the limit on (3 by two orders of magnitude. 
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VI. DISCUSSION AND CONCLUSIONS 



We have obtained the exact energy eigenvalues and eigenstates of the harmonic oscillator 
when the coordinate and momentum operators satisfy the modified commutation relations 
Eq. or 

The energy levels, Eqs. ( p3D and (|5^), grow as for large n. The reason for this 
behavior can be understood as follows: The change of variable from p to p in the ID problem 
changes the kinetic term into a tan^ y^P potential which is bounded at p = ±7r/2A/5. 
For higher dimensions, the effective potential is tan^ V P ~^ P'P plus a centrifugal barrier 
cot^ a//? + P'p which introduces a wall at p = TT/2y/P + j3' in addition to the one at p = 0. 
For higher energy eigenstates, the potentials are in essence square wells, leading to the 
dependence of the energy. Indeed, the energy eigenvalues of a spherically symmetric square 



well potential of radius i: 12^(5 + are given approximately by 

+ + (76) 

The parameter 7, introduced in Eqs. (^) and ( pQ]) has no effect on the energy eigenvalues 
and only results in the wave functions acquiring an extra factor of 

+ + (77) 
which cancels the 7 dependence in the weight function of the inner product, Eqs. (||) and 



The original 

{D-\)\n\ ^^^^ 
fold degeneracy of the n-th energy state is broken, leaving only the 

(D + £-!)! + 3)! 



(D -!)!£! (D- !)!(£- 2)! 



(79) 



fold degeneracy for each value of t due to rotational symmetry . This loss of degeneracy 
can be interpreted as the breaking of self-supersymmetry of the harmonic oscillator [18 . 
The natural question arises whether an analogue exists on the level of field theory as a 
potentially new mechanism for supersymmetry breaking. 

Potential constraints on /? that can be placed by measuring the energy levels of an electron 
trapped in a strong magnetic field have been discussed. Even under optimistic assumptions, 
the constraints that can be imposed are weak. 
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In addition to affecting the energy levels of quantum mechanical systems, the modified 
commutation relations, Eqs. (|I|) and (^), may have other far reaching consequences. In 
subsequent papers, we will discuss their effects on the calculation of the cosmological constant 
P^ , and the motion of macrosopic objects |19| . 
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